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Abstract. Elastic scattering of relativistic electrons from the nucleon yields Lorentz invariant form fac¬ 
tors that describe the fundamental distribution of charge and magnetism. The spatial dependence of the 
nucleon’s charge and magnetism is typically interpreted in the Breit reference frame which is related by 
a Lorentz boost from the laboratory frame, where the nucleon is at rest. We construct a toy model to 
estimate how the charge and magnetic radii of the nucleon are modihed between the Breit and lab. frames. 
This has implications for the ratio of the proton electric to magnetic elastic form factors as a function of 
momentum transfer as well as for determinations of the proton charge radius. Predicted corrections based 
on the model are provided for the rms charge radii of the deuteron, the triton, and the helium isotopes. 
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Elastic scattering of relativistic electrons from the nu¬ 
cleon and nuclei yields their fundamental distributions of 
charge and magnetism in terms of Lorentz invariant form 
factors [l|. In the case of the nucleon, modern CW elec¬ 
tron accelerators using polarization observables have en¬ 
abled the determination of the four elastic proton and 
neutron electric and magnetic form factors with preci¬ 
sion over a large range of momentum transfer 0. This 
increased precision has yielded a number of surprises. The 
ratio of electric to magnetic form factors determined us¬ 
ing polarization transfer disagrees with the previous cross 
section data — the former is seen to decrease dramati¬ 
cally with increasing momentum transfer Q. In fact, while 
viewed as a surprise, this behavior was always expected in 
models based on vector meson dominance (see the dis¬ 
cussions in ii and references therein). Furthermore, 
the charge radius of the proton as determined from pre¬ 
cision elastic electron-proton scattering Q disagrees with 
a high precision determination obtained from the Lamb 
shift in muonic hydrogen Bi- This latter discrepancy 
has sparked considerable interest ilS- In this Letter we 
explore how the nucleon’s charge and magnetism can be 
affected by the choice of reference frame. We note that 
the boost corrections have been previously derived for the 
charge distribution [IMl: here we develop a more gen¬ 
eral model that we apply to both charge and magnetic 
form factors at the largest momentum transfers where 
they have been measured. 
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Consider relativistic elastic electron scattering from 
the nucleon. In single photon exchange approximation, the 
unpolarized elastic eN scattering cross section in the lab. 
system (nucleon at rest) can be written 

( 1 ) 


where the Mott cross section is 


acos0/2 
2EeSin^ 9/2 


( 2 ) 


and the recoil factor is free = EefE'^. Here the incident 
electron has energy Eg, the scattered electron has energy 
Eg and the scattering is through angle 0; also a is the 
fine-structure constant. In all expressions, the extreme rel¬ 
ativistic limit is taken, namely the electron mass is ignored 
with respect to its energy. 

The square of the form factor, E^, can be expressed 
in terms of the electric and magnetic Sachs form factors 
and G^" of the proton and neutron, respectively. 
The normalizations at = 0 are G^(0) = G^(0)//Xp = 
G)(^(0)//x„ = I (the neutron charge form factor is irregu¬ 
lar, but may be fixed with respect to the proton case). As 
usual, one has 

Wi = tGI (3) 

W2 = y^[GI + tGI,], (4) 

where r = |Q^|/4 to^ with rriM the nucleon mass. Using 
these the square of the elastic form factor can be written 
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in several equivalent forms: 


IU 2 -k 2IUi tan^ 0/2 

(5) 


(6) 

vlWl + vtWt, 

(7) 


where the relative flux of longitudinally polarized virtual 
photons e is given by 

e~^ = 1 + 2(1 + t) tan^ 0/2 . (8) 


Using the fact that for elastic eN scattering one has 


P = 




1 


1 + T 


we obtain 


(9) 


VL = P = 


{1 + tY 


,^ = ip + tan^0/2 = ^^^, 


-1 


( 10 ) 

( 11 ) 


The longitudinal and transverse responses are related to 
the lUs or to the Sachs form factors via 


Wl = (1 + t)GI ( 12 ) 

Wt = 2tGIi. (13) 

One definition of a “proton charge radius” — the one 
commonly used in analyses of elastic electron-proton scat¬ 
tering — is obtained through the derivative of G^{Q^) 
with respect to the invariant 4-momentum transfer, Q^, 
namely via 

- dG^m 

dQ^ 

Since both G^{Q^) and are invariants, the so-defined 
quantity with dimensions of length is also. We denote 
this type of radius, as the momentum-space proton 

charge radius. It is a convenient quantity to employ when 
inter-comparing data on elastic ep scattering. What it is 
not is the RMS charge radius of the proton. In principle 
the latter would be found by taking the charge distribution 
of the proton in its rest frame, weighting by r^, integrating 
and taking the square root to obtain the coordinate-space 
charge radius . In some modeling or approaches such 
as lattice QCD one works in coordinate space and proceeds 
in this manner. Non-relativistically the two radii are the 
same, and were the same logic to be applied to heavy 
nuclei, even relativistically the differences would be neg¬ 
ligible. However, the nucleon (and to a lesser extent few- 
body nuclei) is exceptional: the difference between the two 
definitions is larger than the present uncertainties in the 
data. Purists would say that only the momentum-space 
quantity should be used when discussing elastic ep scat¬ 
tering, since this is intrinsically a process where different 
frames are involved. In the lab frame, for instance, the 



proton starts at rest, but must recoil after absorbing the 
momentum transferred via the exchanged virtual photon; 
and choosing the Breit frame does not make the problem 
go away, for then the proton must be moving in both the 
initial and final states (see below). 

Were only electron scattering to be considered, it is 
hard to argue with this point of view, and it is largely our 
non-relativistic bias that motivates us to attempt to de¬ 
fine a quantity that we hope is some sort of radius. Unfor¬ 
tunately, typical modeling is not done in a fully covariant 
way, but rather is done by attempting to provide a picture 
in a specific frame (for instance, through some approach 
such as lattice QCD or employing models such as the bag 
model) and then Fourier transforming to obtain a form 
factor. For any model that is not boostable this procedure 
incurs the same difficulty as the above for the two types of 
radius. Furthermore, as discussed briefly later where elec¬ 
tronic and muonic atomic hydrogen are mentioned, there 
are other measurements involving the concept of a proton 
charge radius that one wants to use to explore whether 
or not they yield radius values that are compatible with 
those found via electron scattering. An important ques¬ 
tion is: which type of radius do these other measurements 
probe? 

Nevertheless, let us not take the point of view that 
only the invariant momentum-space radius defined above 
is relevant and at least attempt to evaluate the potential 
uncertainty in employing a specific definition of the radius 
(later generalized to all four radii for the nucleon) with 
consequence for comparisons of electron scattering deter¬ 
minations with those from the Lamb shifts in electronic 
and muonic hydrogen. We start by summarizing some of 
the basics in the Breit frame. 

The Breit frame is defined by zero energy transfer of 
the virtual photon and reversal of the 3-momentum of the 
target between the initial and final states. Thus, one has 
ujb = 0 and qb = \/\Q‘^\- Since in that frame (with z- 
axis along the momentum transfer vector) one has the nu¬ 
cleon entering with 3-momentum —pB and leaving with 
3-momentum +pb, one has pb = qB/‘2- Thus, the rela¬ 
tivistic 7 -factor relative to the lab. frame for the nucleon 
in that frame is 

7 = Vl -I- r. (15) 

We note that, no matter what reference frame is adopted, 
the nucleon before and after the scattering must have dif¬ 
ferent momenta. In the lab. frame the nucleon begins at 
rest, but must recoil with the full momentum in the final 
state. In the Breit frame (the “brick-wall” frame as dis¬ 
cussed above) it must enter with momentum —qB/2 and 
leave with momentum -\-qB/2. In any case, Lorentz boosts 
are involved and accordingly it is reasonable to expect 
some consequences from these such as Lorentz contrac¬ 
tions and time dilations. For instance, one might expect a 
moving nucleon to be flattened (as nuclei are in relativis¬ 
tic heavy-ion collisions) with corresponding modification 
of its “charge distribution”. In a covariant description in 
momentum space the spinors and operators are typically 
handled appropriately; however, if one tries to interpret 
the Fourier transforms, problems arise and concepts such 
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as the charge and spin distributions are not simply these 
Fourier transforms, as to a large extent they are for heavier 
systems such as nuclei. One might just work in momen¬ 
tum space and never confront these issues, although there 
is strong sentiment that some interpretation should be at¬ 
tempted in coordinate space and in the following we try 
to do this with a “toy model” that, while certainly not a 
fundamental one, is motivated by our intuition. Note that 
scattering involves these problems of reference frames and 
boosts, whereas the Lamb shift does not have the same 
issue for there one has an electron or muon bound to a 
proton with a single frame to consider. Of course, these 
atomic systems have their own issues (relativistic correc¬ 
tions, recoil effects, etc.) 

We now attempt to formulate a non-relativistic version 
of eN scattering. Assume that one puts a single nucleon 
into the lowest level in a very deep harmonic oscillator 
(HO) potential to avoid any recoil problem. In effect, in 
this “toy model” one can imagine there being a very strong 
trap to hold the nucleon essentially at rest; indeed, a nu¬ 
cleon in a heavy nucleus is confined and not allowed to 
carry the full recoil momentum when electrons are scat¬ 
tered from it, providing an example where this situation 
actually occurs. Then using, for instance, the tables of M 
or the review article of [l^ . one can compute the multi¬ 
pole matrix elements of the CO (Coulomb monopole) and 
Ml (magnetic dipole) elastic scattering operators. As dis¬ 
cussed in the standard literature for electron scattering 
from nuclei, and those papers in particular, the differen¬ 
tial cross section may be written as in eq. [T] and [5] where, 
using the non-relativistic limit for the current operators 
together with lsi /2 harmonic oscillator wave functions, 
one obtains the following: 

(lsi /2 ||Mo^°“'|| lsi/ 2 ) = F, (lsi /2 IIMoll lsi/ 2 ) (16) 
(lsi /2 lsi/ 2 ) = —[^"1 (lsi /2 ||Ai((jx)|| lsi/ 2 ) 

— 2{Fi+F2) ■ (lSi/2 ||^U'?^)II 1 ^ 1 / 2 )] ■ (17) 

Here, following 0 , we use the Dirac and Pauli single¬ 
nucleon form factors, as is conventional for non-relativistic 
treatments of electron scattering from nuclei. Either work¬ 
ing directly with the harmonic oscillator wave functions 
and the explicit forms for the current multipole operators 
(see [ 3 ) or (which is easier) using the tables of [l^ one 
finds for the three required reduced matrix elements 

\/47r (lsi /2 ||Mo(gx)|| lsi/ 2 ) = \/2e ^ (18) 

^^(lsi/2 ||^i(gx)|| lsi/ 2 ) = 0 (19) 

V47r(lsi/2 ||A'((gx)|| lsi/ 2 ) = 2e“^. (20) 

Each is proportional to exp(—y) where y = {bq/2)^ with 
b the HO parameter. However, one must multiply by the 
center-of-mass correction which in the non-relativistic HO 
shell model can be computed: it is a multiplicative factor 
of fcm = exp(-|- 2 //A) = exp(-|-j/) for A = 1 and cancels the 
above factor, leaving only the remaining factors obtained 
by using the above-cited tables. Performing the detailed 


developments one can obtain the form factors in the HO 
shell model and, comparing the covariant expressions in 
eq. and m with these results, one can establish a cor¬ 
respondence between the non-relativistic so-defined quan¬ 
tities and their relativistic counterparts 


Non — relativistic 


Relativistic 


Fl = Fi 


F^ = 


2m N 


{Fl + F2)2 


O (l-|-r)G|; (21) 

O 2tGIi (22) 


where Fi and F 2 are the Dirac and Pauli form factors 
of the nucleon. Note that for the transverse form fac¬ 
tor the left-hand expression involves the square of the 3- 
momentum transfer, whereas the right-hand side involves 
the square of the 4-momentum transfer. Equivalently, to 
“de-relativize” the covariant expressions one should define 
the following non-relativistic expressions: 


G nr 

E 

G nr 
M 


= vtt^g 


<rel 


y/T 


-.G 


•rel 
M 5 


(23) 

(24) 


namely, modified in the opposite ways by the factor 7 
(see eq. HH). The motivation is to go from the covariant 
case where boosts (Lorentz contractions and time dila¬ 
tions) must occur, since, in whatever frame one chooses, 
the nucleon must have non-zero momentum at some point 
in the scattering process, to a frame that is more akin to 
our non-relativistic intuition. 

The usual form factor ratio is defined as follows 




G^e 


G\i/ Mp 


(25) 


and is shown in Fig. [T] see Q for references to the data 
and to the so-called GKex vector meson based model @ 
shown as a solid line in the figure. The “de-relativized” 
ratio is then immediately given by 


[Rr = {^ + r)[R,r^ (26) 

and is shown in Fig.[2l Note that this has the boost factor 
squared (going as 1 -I- r and shown in the right panel as a 
red line) rather than just linearly as in the individual form 
factors. Clearly this introduces large modifications at high 
momentum transfers. Indeed, the “de-relativized” results 
are relatively flat as functions of and differ from unity 
by less than roughly 20 %. 

Next, one might attempt to “de-relativize” the defi¬ 
nition of the RMS radii using the same procedures. The 
usual definition is obtained by expanding the jo spherical 
Bessel function for low momentum transfer to obtain 



(27) 


this is simply a re-writing of eq. 1141 That is, we iden¬ 
tify with what was called at the beginning of 
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Fig. 2. (Color online) Left panel: = (1 + r) vs. with data and GKex curve as in Eig. [T] together with 1 + r 

curve referenced to the right-hand axis. 
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the paper. We then proceed to “de-relativize” to obtain a 
quantity that might be more intuitive, based on our toy 
model. For this one must include the \/l + r factor in 
eq. 1231 obtaining at small momentum transfer 


Gl^ ={1 + -T 


■ [teT + 


= 1 - 


2 2 r re/1 ^ ^ 

[te] - - 


r + • ■ 


= 1 - -Tm% + • 

and leading to the relationship 


(28) 

(29) 

(30) 


r 


nr _ 

E — 



(31) 


where one has 

Ap = = 0.0332 W (32) 

4mp 

An = = 0.0331 fm^ (33) 


for protons and neutrons, respectively. The same argu¬ 
ments for the magnetic form factor where the required 
boost factor is now l/-\/l + t (see eq. [24]) leads to the 
expression 

+ A (34) 

We note that the same result for the proton charge form 
factor was obtained in jl4 , arguing from a very different 
point of view: see also [ll| . on which that work is based. 

It is important to understand that these differences be¬ 
tween relativistic and non-relativistic radii do not go away 
if electron scattering data are obtained at ever smaller 
values of the momentum transfer. As the above expres¬ 
sions clearly show, the relativistic boost factor arising from 
(1 -|-t)^^/^ deviates from unity at order however, that 
is the order needed to extract the charge or magnetic radii. 
In other words, being locked together at the same order 
when expanding in powers of the effects can never be 
separated, no matter how small the momentum transfer 
becomes. 

Specifically, usiig the Bernauer value for the proton 
rms charge radius and the PDG values [l^ for the 
proton magnetic, neutron charge, and neutron magnetic 
rms radii, one has the following: 


= ^7(0.879)2 - 0.0332 
= 0.860 ±0.008 fm 
r^^p = ^(0.777)2 + 0.0332 

= 0.7981 ±0.013 ±0.010 fm 
= v'(0.862)2 -h 0.0331 

= 0 . 8810 l°;r 8 fm 

= -0.1161- 0.0331 
= -0.1492 ± 0.0022 W . 


(35) 

(36) 

(37) 

(38) 


Here the uncertainties are taken from the Bernauer work Q 
and from the PDG compilation [l^, respectively. Note 
that the electric result for the neutron is traditionally ex¬ 
pressed as the square of the radius, which is negative. 

What has been obtained above from electron scatter¬ 
ing through the simple toy model are radii that start from 
the quantities measured in electron scattering (r™"™ = 
r’'®^) to quantities that are more like model 

employed in this work thus leads us to make the identifi¬ 
cation o r"’’. Namely, when coordinate-space radii 

are wanted, rather than the momentum-space ones mea¬ 
sured in electron scattering, the toy model motivates tak¬ 
ing into account the corrections discussed above. While 
the model is not “fundamental”, at least the fact that the 
relativistic and non-relativistic quantities so-defined are 
different should be cause for some concern that the con¬ 
cept of a charge or magnetic radius is not totally robust. 

Let us take this one step further and bring in the 
values of the “proton charge radius” determined via the 
Lamb shifts in electronic and muonic hydrogen which are 
known to have some charge-distribution dependence be¬ 
cause the proton is not a point particle, but has a finite 
extent. Note that these atomic systems are described in 
coordinate space and not, as for electron scattering, in 
momentum space. In the course of developing the formal¬ 
ism for studies of atomic hydrogen it is natural to invoke 
the Fourier transform of the proton charge distribution to 
inter-relate what is desired in this case with what is mea¬ 
sured in electron scattering; however, this is the essence 
of the issue. Namely, it would appear to be more natu¬ 
ral that the Lamb shift problem involves the coordinate- 
space radius and not its momentum-space analog, 

i.e., to use r'^p-, not its relativistic partner. One might be 
tempted to use the fact that the effective value of QrE,p 
is very small for these atomic systems (~ 10“^); however, 
the above argument on the “locking” of the boost factor 
with the radius shows that the smallness of this product 
is not sufhcient to make relativistic and non-relativistic 
radii effectively the same. It is not the scale of momentum 
transfer that is critical (as long as it is small enough to 
allow only terms of quadratic order to be considered), but 
the fact that a scattering process and measurements of an 
atomic system are more naturally studied in momentum 
space and coordinate space, respectively. We assume that 
analysis of the Lamb shift entails using the coordinate- 
space version of the proton charge radius which, in our 
toy model, means r'^p. 

The proton charge radius discrepancy has arisen from 
the different values resulting from a precise determina¬ 
tion using the Lamb shift in muonic hydrogen (0.84087 
± 0.00026 ± 0.00029) [1] disagreeing with the GODATA 
2010 value (0.8775 ± 0.0051) [2^, largely determined by 
the most precise value resulting from elastic electron pro¬ 
ton scattering This amounts to more than 4% differ¬ 
ence, whereas the stated total uncertainty in the Bernauer 
electron scattering result is quoted as 0.9%. On the other 
hand, the correction resulting from the boost between 
Breit and lab. frames as calculated in our toy model de¬ 
creases the electron scattering determination of rE,p to- 
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Fig. 3. Timeline of recent determinations [61. IflI. I2,T12 h| of the proton charge radius rE,p- For earlier experimental results see [7|]. 
The Bernauer value for the proton charge radius, namely is indicated with a black cross lying at 0.879, while 

its corrected value is indicated by the white cross at 0.860. The shaded bands show the Icr and 2a uncertainties 

about the latter and may then be compared with the Lamb shift value indicated by a star and lying at about 0.84. 


wards the muonic hydrogen value. The resulting discrep¬ 
ancy in the different determinations of rE,p is cut in half 
using the corrected value and now differs from the muonic 
Lamb shift value by only about 2%. Further, the corrected 
value here for the proton charge radius is not inconsis¬ 
tent with the value determined using the hydrogen atom, 
within experimental uncertainty. Fig. |2| shows a selection 
of determinations of the proton rms charge radius since 
the year 2000 including all the work cited here. In partic¬ 
ular, it shows the electron scattering determination from 
Bernauer and our corrected value with both Icr and 2cr 
uncertainties. If one accepts the ideas we have developed 
in this paper, then one is left with only a small amount 
to explain and, the remaining discrepancy being commen¬ 
surate with the size of a, might lead one to question the 
other model dependences in the problem (experimental 
systematic uncertainties, radiative corrections, etc.) 

The independent recoil correction to the rms charge 
radius derived from elastic electron scattering can also be 
evaluated for the deuteron, the triton, and the helium iso¬ 
topes. For heavier nuclei, the correction factor above is 
replaced as follows: 


3 . 

4:171% \mAj ' 


(39) 


where is the mass of that heavier system. In Table [TJ 
we summarize the results based on our model. We note 
that the corrected deuteron rms charge radius from elec¬ 
tron scattering is in excellent agreement with the recent 
high precision value obtained from measurements of the 
Lamb shift in muonic deuterium. We provide predictions 
for the rms charge radii of the helium isotopes, which are 
being determined to high precision in ongoing experiments 
that employ measurement of the Lamb shift in muonic 
atoms. New high precision measurements of the Lamb 
shift in electronic hydrogen are in progress and results 
are expected soon. 

We have also considered how the conclusions here can 
be validated by experiment. We point out that the correc¬ 
tion we derive cannot be separated by going to lower in 
electron scattering experiments nor by comparison of elas¬ 
tic electron and muon scattering on the proton. However, 
based on the estimates made here, a precision comparison 
of electron scattering from the proton with electron scat¬ 
tering from the deuteron at low should deviate at the 
level of about 1.5%, since the corrections we derive dif¬ 
fer at this level. For highest precision, such measurements 
should be carried out using the same apparatus and sys- 
tematics must be minimized. We note that internal radia¬ 
tive corrections, which arise mainly from the incident and 














T.W. Donnelly et al.: Interpretations of Elastic Electron Scattering 


7 


Nucleus 

te 

e-scattering 

(fm) 

ve 

e-scattering corrected 
(fm) 

ve 

muonic atom 
(fm) 


0.879 ± 0.008 

0.860 ± 0.008 

0.84087 ± 0.00039 


2.13 ± 0.01 [3^ 

2.12 ± 0.01 

2.12562 ± 0.00078 [31J 


1.755 ± 0.087 [3^ 

1.751 ± 0.087 


®He 

1.959 ± 0.034 [3^ 

1.955 ± 0.034 


^He 

1.680 ± 0.005 [3^ 

1.678 ± 0.005 



Table 1. Comparison of charge radii measured by electron scattering from various light nuclei and the corrected values following 
the procedure proposed herein. The muonic determinations for the proton and the deuteron (preliminary value - see M) are 
also given. We note that high precision values for the helium isotopes will be forthcoming from muonic atom experiments. 


scattered electrons, should be quite similar for proton and 
deuteron targets. 

In summary, we argue that corrections between the 
Breit and lab. frames are important in interpreting the 
form factors of the nucleon as determined in relativistic 
electron scattering. We have constructed a toy model to 
estimate these corrections. In this toy model, the observed 
significant decrease of the ratio of the proton elastic form 
factors as a function of is understood as a predomi¬ 
nantly relativistic effect. Furthermore, in this toy model, 
the proton charge radius as determined in electron scat¬ 
tering has a correction that reduces its value towards that 
resulting from the precision determination using the Lamb 
shift in muonic hydrogen. The toy model provides predic¬ 
tions for the rms charge radii of the deuteron, triton, and 
helium isotopes. 

The ideas presented here should not be viewed as pro¬ 
viding a detailed model of how boost effects might enter 
in electron scattering from the nucleon, but as a warning 
that they could play a role in electron scattering, but not 
in determinations of the Lamb shift. Much more satisfy¬ 
ing would be to have a model of the nucleon that could 
be boosted, although it is not obvious that such a model 
exists at present. For instance, in models where relativis¬ 
tic quarks are confined via a “bag” one must confront the 
problem of how to boost the latter in going between the 
frames that inevitably enter in electron scattering. One 
possible, but different, related case that could be studied 
to test the ideas is that of relativistic (covariant, boost- 
able) modeling of the deuteron [s^- There, on the one 
hand, one could directly compute the elastic form factor, 
while on the other hand one could compute the ground- 
state charge distribution and then Fourier transform it 
to momentum space. Upon comparing the two results it 
is likely that differences will be found that relate to the 
boost issues raised in the present study. 

We thank Jan Bernauer for valuable discussion. The authors’ 
research is supported by the Office of Nuclear Physics of the 
U.S. Department of Energy under grant Contract Numbers 
DE-SC0011090 and DE-FG02-94ER40818. 
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